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Abstract: 

We examine the contribution of Multiple Parton Interactions to Z + n~jets production at 
the LHC, n = 2,3, 4, where the Z boson is assumed to decay leptonically. 
We compare the results obtained with the correlated GS09 double parton distribution 
function with those obtained with two instances of fully factorized single parton distribution 
functions: MSTW2008LO and CTEQ6L1. 

It appears quite feasible to measure the MPI contribution to Z+2/3/4 jets already in the 
first phase of the LHC with a total luminosity of one inverse femtobarn at 7 TeV. If as 
expected the trigger threshold for single photons is around 80 GeV, Z + 2-jets production 
may well turn out to be more easily observable than the 7 + 3-jets channel. The MPI 
cross section is dominated by relatively soft events with two jets balancing in transverse 
momentum. 
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1. Introduction 

The Q CD-improved Parton Model forms the basis of our understanding of high-energy 
hadron scattering. In this framework each hadron is described as a collection of essentially 
free elementary constituents. The interactions between constituents belonging to different 
colliding hadrons are the seeds of the complicated process which eventually leads to the 
particles observed in the detector. In this conceptual scheme it is quite natural to envisage 
the possibility that more than one pair of partons might interact in a single hadronic 
impact. This kind of events are referred to as Multiple Parton Interactions (MPI) while 
those in which only a single pair of partons produce a hard scattering are described as 
Single Parton Interactions (SPI). 

Multiple Parton Interactions have been detected in high transverse momentum hadron 
collisions both at the ISR at CERN [1] and at the Tevatron at Fermilab [2, 3, 4]. The 
measured cross sections imply that MPI could provide a non-ncgiigiblc background to all 
sort of interesting reactions since MPI rates at the LHC are expected to be large. At 
smaller transverse momentum MPI have been shown to be necessary for the successful 
description of the underlying event both in Pythia [5, 6, 7] and in Herwig [8, 9]. The wide 
range of phenomena in which MPI are involved highlights the urgency of a more thorough 
understanding of these reactions both experimentally and from a theoretical point of view. 

The theoretical investigation of MPI has a long history [10, 11, 12, 13, 14] and have 
experienced a renewed interest in more recent times [15, 16, 17, 18, 19, 20, 21, 22, 23]. 

The basic formalism can be readily described starting from the standard expression 
for the SPI cross section: 



where ti is the factorization scale which characterizes the interaction and at which the 
parton distribution functions Fi(xi,ti) for parton i to have momentum fraction xi are 




(1.1) 
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evaluated. Eq.(l.l) can be rewritten in term of parton distributions which depend on the 
transverse coordinates as well as on the longitudinal momentum fraction as: 

(T^A) =Y1 I '^iixi,bi,h)afj{xi,yi)rj{yi,bi- P,ti)dxidyidHid'^P (1.2) 

where (5 is the usual impact parameter. Making the reasonable assumption that the de- 
pendence on the momentum fraction and that on the transverse position factorize 

Ti{x,h) = E,{x)x f{h) (1.3) 

and that the latter is a universal function for all kind of partons fixes the normalization of 
the transverse distribution: 

j f{h)f{h-l3)d%d^P = j T{P)d^p = l (1.4) 

where we have defined the overlap function T(/3) = J f{b)f{b — (3) cPb. 

Analogously we can write the Double Parton Interaction (DPI) cross section as follows: 

'^iA,B) = / ^i,ji^i^bi,ti,X2,b2,t2)a^kixi,yi)crfA^2,y2) (1-5) 

''^''''^ xrk,iiyi,bi- /3,ti,y2,b2- P,t2)dxidyidHidx2dy2d%d'^P 

where ti and t2 are the factorization scales of the two scatterings; m is a symmetry factor 
which is equal to one if the reactions A and B are identical and equal to two if they are 
not. 

Separating the transverse part, rij{xi,bi,ti,X2,b2,t2) = Fi,j{xi,ti,X2,t2) x /(6i) x 
7(62) Eq.(1.5) becomes 

'^{A,B) = 21 7 J Fi,jixi,ti,X2,t2)a^kixi,yi)crf,i{x2,y2) (1-6) 

''J'*^'' X Fk^i{yi,ti,y2, t2) dxi dyidx2 dy2 

where ^ ^. 

' T\l3)d^l3. (1.7) 



If one makes the further assumptions that double parton distributions reduce to the 
product of two independent one parton distributions, Fjj = x Fj, the DPI cross section 
can be expressed in the simple form 

This last assumption however, even though rather common in the literature and quite 
convenient from a computational point of view, is clearly incorrect. In Ref. [15, 16] it was 
shown that correlations between the value of the double distribution functions for different 
values of the two momentum fractions xi,X2 are to be expected, even under the assumption 
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of no correlation at some scale Qq, as a consequence of the evolution of the distribution 
functions to a different scale Q, which is determined by an equation analogous to the usual 
DGLAP equation [24, 25, 26]. 

A large number of studies have evaluated the MPI contribution to several high en- 
ergy processes [27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40], including Higgs and 
electrowcak vector boson production. Other studies [41, 42, 43] have in addition focused 
on the differences between final states produced in SPI and in MPI as a tool to reduce 
the background due to MPI or alternatively to separate MPI processes from SPI ones 
and gain more detailed experimental information on Multiple Parton Interactions. All the 
aforementioned studies have assumed complete factorization of double Parton Distribution 
Functions (dPDF). 

In [44] the corrections to the factorized form for the dPDF have been estimated. They 
depend on the factorization scale, being larger at larger scales Q, and on the x range, 
again being more important at larger momentum fractions. For Q = Mw and x ^ 0.1 the 
corrections are about 35% for the gluon-gluon case. 

Recently Gaunt and Stirling [45] have developed a set of dPDF which satisfy a collec- 
tion of momentum and number sum rules. These sum rules are preserved by the evolution 
equations [15, 16] and therefore are obeyed at any scale Q once they are satisfied at an input 
scale Qq. The GS09 set is based on the MSTW2008LO single Parton Distribution Func- 
tions (sPDF) [46]. Gaunt and Stirling also provide a program which evolves the dPDF from 
the input scale to any scale and a set of publicly available dPDF grids. In their published 
form the GS09 set deals with the case of two identical scales ti and t2 in the distribution 
functions Fij, but this limitation has been recently dropped. 

In Rcf. [47] the GS09 dPDF have been employed in a study of same-sign W pair 
production at the LHC, including the background due to W^Z{'y*) production. From 
the ratio R = Aay^+w+^w-w- / '^w+w- ' "'^l^ich is equal to one for factorized dPDF, a 
violation of factorization at the 20% to 30% level is reported. 

In this paper we examine the contribution of MPI to Z+n-]e\,s production at the LHC, 
n = 2, 3, 4, where the Z boson is assumed to decay leptonically. These processes have the 
advantage of a much larger cross section than same-sign WW production and therefore are 
more likely to allow detailed studies of MPI at the low luminosity, about 1 fb~^, foreseen 
for the first two years of operation at the LHC with ^/s = 7 TeV. While the cross section 
for Z + nj is smaller than for W + nj because of the smaller leptonic branching ratio, 
the former is cleaner from an experimental point of view since isolated, high pT charged 
leptons, which are the main tool for W detection, can be copiously produced in B-hadron 
decays [48, 49, 50] while no comparable mechanism exists for generating lepton pairs of 
mass in the Mz region. As pointed out in [47] Z{'y*) + jets production, with one of the 
leptons undetected, can also mimic W + nj processes. 

Z + nj production probes initial state parton combinations which are different from 
those probed in 1^=*=!^='= processes. The latter, at lowest order, are always initiated by 
four-fermion states, mainly udud. The former, on the contrary, typically have at least two 
gluons in the initial state since the largest component [42, 43] involves a two jet process 
which is dominated by gluon-gluon scattering. 
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For comparison we also present the predictions for 7 + 3j production, the reaction from 
which the most recent and precise estimates of Ceff have been extracted at the Tevatron. 
This measurement will undoubtedly be performed again at the LHC [51]. 

NLO QCD corrections are or will soon be available for all SPI processes leading to an 
electroweak vector boson in association with up to four jets [52, 53, 54, 55, 56, 57, 58, 59]. 
The Drell-Yan cross section is known at NNLO [60]. Measurements at the Tevatron show 
good agreement between NLO calculations and data [61, 62]. These new developments 
open the possibility of validating the predictions using events with large visible energy, 
where the MPI contribution is small, and then using them for a direct measurement of 
the MPI cross section at smaller total invariant masses in parallel with more data driven 
analysis similar to those of CDF and DO. 

In the following we compare the results obtained with the GS09 dPDF with those 
obtained with two instances of fully factorized sPDF: MSTW2008LO [46] and CTEQ6L1 
[63]. Hence we can estimate, even in the absence of a proper dPDF set based on CTEQ6, 
the dependence of MPI predictions on the choice of PDF, a study that to our knowledge 
has not been performed before. 

We have considered three center of mass energies for the LHC: ^/s = 7 TeV, \/s = 10 
TeV and -y/s = 14 TeV. This allows us to study the properties of MPI processes while the 
relevant range of momentum fractions for the dPDF shifts to smaller values as the energy 
increases. 

Given the strong similarities between the production mechanism of Z + jets, W + jets 
and 7 + jets we expect that the conclusions reached in the present paper for Z + jets 
production concerning the ratio of the MPI to the SPI contribution, the effect of correlations 
in MPI and the dependence on the PDF choice will be applicable also io W + jets and 
7 + jets production. 

We will confine ourselves to Double Parton Interactions and neglect Triple and Higher 
Order Parton Interactions. Triple Parton Interactions are expected to be significantly less 
abundant than Double ones, even though it has been argued that they could be indeed 
detected at the LHC [42, 43]. 

In Sect. 2 the main features of the calculation are discussed. Then we present our 
results in Sect. 3. Finally we summarize the main points of our discussion. 

2. Calculation 

The MPI processes which contribute at leading order to Z + n-jets through Double Parton 
Interactions are those in which an event producing k jets is superimposed to an event 
producing a Z-boson and (n — k) jets, k = 2, . . . ,n. 

At the Tevatron, CDF [2, 3] has measured erg// = 14.5± 1.7l2'3 a value confirmed 
by DO which quotes a^ff = 15.1 it 1.9 mb [4]. In Ref. [17] it is argued, on the basis of 
the simplest two channel cikonal model for the proton-proton cross section, that a more 
appropriate value at ^/s = 1.8 TeV is 10 mb which translates at the LHC into (^^ff^ ~ ^"^ 
mb. Treleani then estimates the effect of the removal by CDF of TPI events from their 
sample and concludes that the CDF measurement yields Cg// ~ 11 mb at Tevatron energies. 
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In the following we use cr^ff = 12.0 mb for all LHC center of mass energies, with the 
understanding that this value is affected by an experimental uncertainty of about 15% 
and that it agrees only within 30% with the predictions of the eikonal model. Since cie// 
appears as an overall factor in our results it is easy to take into account a different value. 

It is worth mentioning that at present there is a discrepancy between the value of 
(Te/ / extracted by CDF and DO and the one which is effectively employed by Pythia whose 
normalization is derived mainly from comparisons with small p-p data which dominate the 
total cross section. The description of MPI in PYTHIAS [64] assumes that interactions can 
occur at different pt values independently of each other inside inelastic non-diffractive 
events. The expression for a DPI cross section becomes therefore: 

O" =< fimpact > <7l • CF2/(TND/k (2.1) 

where fjjv/) is the total non-diffractive cross section and fimpact is an enhancement/depletion 
factor chosen event-by-event to account for correlations introduced by the centrality of the 
collision. This quantity is typically averaged during an entire run to calculate < fimpact > 
in Eq. 2.1. Typical values at the center of mass energy of 10 TeV are 1.33 for < fimpact > 
and 51.6 mb for a^D- Comparing Eq. 2.1 with Eq. 1.8 tells us that PYTHIAS predicts an 
effective (Jeff=crND/< fimpact > which is about a factor three larger than the one actually 
measured at the Tevatron. I believe that this issue deserves careful consideration and that 
new measurements of high py MPI reactions would be quite welcome. 
All samples have been generated with the following cuts: 



PT, > 30 GeV, \7]j\< 5.0, 

PT, > 20 GeV, \r,e\ <2.5, (2.2) 
PT, > 30 GeV, \r]^\ < 2.5, 
ARjj > 0.1 , ARji > 0.1 , ARj^ > 0.1 

where j = u, u, d, d, s, s, c, c, b, b, g and / = e+, e~ , fi~ , 

The Z + 4-jets sample has been generated with PHANTOM [65, 66, 67], while all other 
samples have been produced with MADEVENT [68, 69]. Both programs generate equal weight 
events in the Les Houches Accord File Format [70] . All samples have been generated using 
CTEQ6L1 [71] parton distribution functions. The QCD scale (both in ag and in the parton 
distribution functions) has been taken as 

n 
i=l 

where n is the number of final state partons, for all reactions with the exception of qq 
Z+Z" for which the scale has been set at = M^. The scale in Eq.(2.3) is similar, though 
not identical, to the scale advocated in Refs. [58, 59] for vector boson plus jets production 
at NLO. 

The results shown in the following under the CTEQ heading have been obtained com- 
bining at random one event from each of the reactions which together produce the desired 
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final state through MPI. When needed, we have required that each pair of colored par- 
tons in the final state have a minimum AR separation. This implies that the combined 
cross section does not in general correspond to the product of the separate cross sections 
divided by a^ff because the requirement of a minimum separation for all jet pairs induces 
a reduction of the cross section when additional pairs are formed in superimposing events. 

The results shown under the MSTW and GS09 headings have been obtained through 
a reweighting procedure by the appropriate ratio of parton distribution functions and 
coupling constants. For instance, an event like {qiQi — )• gl~^l~) {gg — > gg), constructed 
from two events generated separately with CTEQ6 PDF, can be transformed in a weighted 
event with MSTW2008 PDF multiplying its original weight by 



MSTW, . MSTW, s MSTW,, , j^MSTW, . MSTW, , MSTW,, ,9 

j^_ Fi ih)F. ih) (ti) {t2)Fg (^2) a, {12? 

where ti, t2 are the factorization scales for q^qi — )• gl'^l~ and gg — )■ gg respectively. The 
factorization scales have been read off from the event files. The second and fourth factors 
in Eq.(2.4) take into account the different values of the strong coupling constants for the 
two different sets of PDF: a^Jo'^(Mz) = 0.130 while a"lo {Mz) = a°^°\Mz) = 0.139. 
The only difference for the GS09 case would be that the correlated dPDF Fij(ti,t2) would 
appear instead of the uncorrelatcd product Fi{ti)Fj{t2) and so on. The resulting events 
are no longer unweighted. The error on the cross section introduced by the reweighting 
procedure is essentially negligible because of the large size, about 5 x 10^ events, of the 
samples. Reweighting can also be employed to estimate the sensitivity of our tree level 
results to variations of the scale Eq.(2.3) using a straightforward modification of Eq.(2.4). 

All results are obtained with the following values for the clcctrowcak input parameters: 
Mz = 91.188 GeV, Mw = 80.40 GeV, Gp = 0.116639 x 10"^ GeV'^. 

We work at parton level with no showering and hadronization. Color correlations 
between the two scatterings have been ignored. 



3. Results 

The total cross sections for SPI and DPI production for Z-|-2-jets, Z+S-jets and Z-|-4-jets 
are presented in Tab. 1, Tab. 2 and Tab. 3 respectively. In all cases the cuts in Eq.(2.2) have 
been imposed. The analysis has been repeated requiring a larger separation between jets; 
the results for ARjj = 0.5 and for ARjj = 0.7 arc also shown. In our estimates below we 
have only taken into account the muon decay of the Z boson. The Z — t- e^e~ channel gives 
the same result. The possibility of detecting high pT taus has been extensively studied in 
connection with the discovery of a light Higgs in Vector Boson Fusion in the t'^t~ channel 
with extremely encouraging results [72]. Therefore, the tau decay of the Z can be expected 
to contribute significantly to MPI studies. 

The total cross sections for SPI and DPI production for 7 + 3-jcts are shown in Tab. 4 
with increasing jet-jet separation. It should be mentioned however that at the LHC trigger 
thresholds for single photons are foreseen to be much higher than those for double leptons 
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14 TeV 


10 TeV 


7 TeV 


Z + 2j 


CTEQ 


MSTW 


GS09 


CTEQ 


MSTW 


GS09 


CTEQ 


MSTW 


GS09 


ARjj = 0.1 


SPI 


56.71 


65.35 




33.11 


37.97 




17.97 


20.48 




DPI 


11.27 


11.37 


15.50 


1.80 


6.35 


G.G8 


1.88 


2.01 


2.66 


ARjj = 0.5 


SPI 


52.65 


60.70 




30.63 


35.15 




16.56 


18.88 




DPI 


11.27 


14.37 


15.50 


4.80 


6.35 


6.68 


1.88 


2.61 


2.66 


ARjj = 0.7 


SPI 


51.53 


59.41 




29.95 


34.38 




16.17 


18.45 




DPI 


11.27 


14.37 


15.50 


4.80 


6.35 


6.68 


1.88 


2.61 


2.66 



Table 1: Z + 2 jets, Z fi^fi cross sections in pb. Cuts as in Eq.(2.2) with increasing angular 
separation between jets, ^Rjj- 





14 TeV 


10 TeV 


7 TeV 


Z + 3j 


C^TEQ 


MSTW 


GS09 


C^TEQ 


MSTA\' 


GS09 


CTEQ 


MSTA\' 


GS09 


^Rjj = 0-1 


SPI 


21.62 


26.25 




11.75 


14.18 




5.77 


6.91 




DPI 


2.93 


4.06 


4.20 


1.10 


1.61 


1.60 


0.37 


0.58 


0.55 


ARjj = 0.5 


SPI 


15.71 


19.10 




8.46 


10.23 




4.11 


4.93 




DPI 


2.70 


3.75 


3.88 


1.02 


1.49 


1.48 


0.34 


0.54 


0.51 


ARjj = 0.7 


SPI 


14.13 


17.18 




7.59 


9.18 




3.67 


4.41 




DPI 


2.59 


3.60 


3.73 


0.97 


1.43 


1.42 


0.33 


0.52 


0.49 



Table 2: Z + 3-jets, Z cross sections in pb. Cuts as in Eq.(2.2) with increasing angular 

separation between jets, ARjj. 



[73, 74, 75]. While pair of leptons are expected to be triggered on for transverse momenta of 
about 15 GeV, single photons will be detected only when their transverse momenta is larger 
than about 80 GeV at the design energy of 14 TeV. At lower energies and instantaneous 
luminosities the threshold could be smaller. Even at design luminosity and center of mass 
energy a lower threshold could be allowed with some pre-scaling. Since MPI processes 
are known to decrease sharply with increasing transverse momenta, we present in Tab. 5 
the predictions for pr^ > 80 GeV while the results in Tab. 4 are mainly intended for low 
luminosity data taking. 

The Single Particle Interaction MSTW results are larger than those obtained with the 
CTEQ PDF by an amount which varies between 15% for Z+2j to 27% for Z+Aj, increasing 
as expected with the power of ag in the amplitude. The Double Particle Interaction MSTW 
results are larger than those obtained with the CTEQ PDF by an amount which varies 
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14 TeV 


10 TeV 


7 TeV 


Z + 4j 


CTEQ 


MSTW 


GS09 


CTEQ 


MSTW 


GS09 


CTEQ 


MSTW 


GS09 


ARjj = 0.1 


SPI 


8.80 


11.16 




4.23 


5.33 




1.80 


2.25 




DPI 


1.21 


1.92 


1.82 


0.41 


0.71 


0.G6 


0.12 


0.23 


0.20 


ARjj = 0.5 


SPI 


4.26 


5.41 




2.00 


2.53 




0.83 


1.04 




DPI 


0.96 


1.53 


1.50 


0.33 


0.56 


0.52 


0.10 


0.18 


0.16 


ARjj = 0.7 


SPI 


3.34 


4.24 




1.56 


1.97 




0.64 


0.80 




DPI 


0.87 


1.39 


1.35 


0.29 


0.51 


0.47 


0.09 


0.16 


0.14 



Table 3: Z + 4-jcts, Z fi^fi cross sections in pb. Cuts as in Eq.(2.2) with increasing angular 
separation between jets, ^Rjj- 



between 30% and 90%. The larger shift is due to the smaller scales for the two individual 
scatterings compared to a single interaction event with the same final state particles. The 
predictions for the GS09 correlated dPDF are larger than those with MSTW uncorrelated 
ones for -y/s = 14 TeV and y/s = 10 TeV while they are smaller for -y/i = 7 TeV. The 
difference is at most of 15%. Taking into account the errors in the measurement of fg// 
we conclude that the uncertainties due to the choice of PDF and to correlation effects are 
reasonably under control. 

These variations should be compared with the uncertainty due to scale variation in PDF 
and in the strong coupling constant. In order to estimate the latter we have reweighted our 
samples changing the scale in Eq.(2.3) by a factor of two in either direction for two limiting 
cases, namely Z + 2j production at ^/s = 7 TeV and Z + 4j production at -^/s = 14 TeV. 
In both instances we have used MSTW PDF and ARjj = 0.5. For Z + 2j production at 
^/s = 7 TeV the cross section changes by +14%/- 13% when the scale is halved /doubled; 
in the case of Z + Aj production at ^/s = 14 TeV the corresponding shifts are +57%/-29%. 
The processes we are interested in therefore are not overly sensitive to scale variations. 
The corresponding uncertainty is of the same order than that related to PDF choice. 

The effects of higher order corrections are more difficult to estimate since no NLO 
calculation for MPI processes is available. QCD one loop calculations arc available for 
vector boson production with up to four jets [52, 53, 54, 55, 56, 57, 58, 59] and arc typically 
of order 10% with the exception of Drell-Yan inclusive production [76] where they are of 
the order of 50%. NLO corrections for the inclusive jet cross section at the LHC have been 
presented in Ref . [77] . For small transverse momenta, as the ones we are interested in this 
paper, they are of the order of 10%. 

The ratio between the MPI and SPI cross sections increases with the collider energy, 
that is with decreasing average momentum fractions carried by the incoming partons. It 
also increases with the ARjj separation because of the absence of correlations between the 
final state partons originating in the independent scatterings which compose MPI events. 
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14 TeV 


10 TeV 


7 TeV 


7 + 3j 


CTEQ 


MSTW 


GS09 


CTEQ 


MSTW 


GS09 


CTEQ 


MSTW 


GS09 


ARjj = 0.1 


SPI 


5921.7 


7341.4 




3484.2 


4302.5 




1884.9 


2317.1 




DPI 


436.9 


612.7 


663.2 


176.4 


262.4 


273.8 


61.4 


103.4 


103.0 


ARjj = 0.5 


SPI 


4516.7 


5610.2 




2637.2 


3263.8 




1415.8 


1744.6 




DPI 


422.2 


593.7 


642.2 


170.3 


254.1 


264.9 


62.0 


100.0 


99.4 


ARjj = 0.7 


SPI 


4137.4 


5142.9 




2411.1 


2986.4 




1290.3 


1591.3 




DPI 


407.5 


574.6 


621.0 


164.1 


245.7 


255.9 


59.7 


96.6 


96.0 



Table 4: 7 + 3-jets cross sections in pb. Cuts as in Eq.(2.2) with increasing angular separation 
between jets, ARjj. 





14 TeV 


10 TeV 


7 TeV 


7 + 3j 


CTEQ 


MSTW 


GS09 


CTEQ 


MSTW 


GS09 


CTEQ 


MSTW 


GS09 


Ai?,, = 0.1 


SPI 


944.8 


1142.4 




524.0 


629.9 




256.8 


306.6 




DPI 


18.4 


29.0 


29.0 


6.84 


11.88 


11.22 


2.17 


4.24 


3.72 


ARjj = 0.5 


SPI 


671.5 


813.0 




368.30 


443.38 




177.4 


212.0 




DPI 


17.7 


28.1 


28.1 


6.59 


11.49 


10.85 


2.09 


4.09 


3.58 


ARjj = 0.7 


SPI 


599.7 


726.4 




328.18 


395.29 




156.6 


187.2 




DPI 


17.1 


27.2 


27.2 


6.34 


11.10 


10.47 


2.00 


3.94 


3.45 



Table 5: 7 + 3-jets cross sections in pb. Cuts as in Eq.(2.2) and pr^ > 80 GeV, with increasing 
angular separation between jets, ARjj. 



For Z + nj processes and taking ARjj = 0.5 as an example, the ratio is of the order of 
10% for ^ = 7 TeV and grows to about 25% at ^/s = U TeV. The results for 7 + 3-jets 
show a similar behaviour with somewhat smaller fractions of MPI events to SPI ones which 
however depend drastically on the px^ cut. For > 30 GeV they range between 5 and 
10% while for px^ > 80 GcV they are at the percent level. 

If we consider the MPI processes as our signal and the SPI ones as the corresponding 
background, we can estimate the prospect of measuring MPI in a given final state from the 
standard S/s/B significance. Using for S the result obtained with GS09 PDF and for B 
the result for the MSTW set and assuming a luminosity of one inverse femtobarn at 7 TeV, 
the significancies extracted from Tabs. 1-3, in the Z — t- n'^fi^ channel alone, are 19/7/5 
for Z+2/3/4 jets with ARjj = 0.7. The corresponding number of expected MPI events 
are 2600/500/140. Therefore it appears quite feasible to measure the MPI contribution to 
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Figure 1: On the left: distribution of tlie angular separation in the transverse plane between the 
two highest px jets in Z + Aj events. On the right: distribution of the total visible mass, pO^j 
in ^ + 2j events. For both plots ^/s = 7 TeV, ARjj = 0.7. 



Z+2/3/4 jets already in the first phase of the LHC. 

The significance of 7 + 3-jets depends on the trigger strategies. If the threshold for sin- 
gle photon detection can be brought in the 30 GeV range then the much larger production 
rate, about ten times that of -^(/U/x) + 2j, provides the best opportunity for an early mea- 
surement of MPI at the LHC. If, on the contrary, the photon trigger cannot substantially 
deviate from about 80 GeV, Z + 2j production looks more promising than the 7 -|- 3-jets 
channel whose significance becomes similar to that of Z + 3j. Anyway, in order to go 
beyond measuring a^ff and start to extract the double parton distribution functions from 
the data, one should measure the MPI fraction of as many channels as possible, exploiting 
the fact that different reactions are initiated by different combinations of partons. 

The contribution to the MPI Z + n-jets cross section due to two jet production in 
association to Z + (n — 2)-jets processes is in all instances the largest one, therefore, even 
with more than two jets in the final state, the majority of MPI events are expected to 
contain a pair of jets which are back to back in the transverse plane. This is confirmed by 
the left hand side of Fig. 1 which displays the distribution of the angular separation Acf) 
between the two highest jets in Z -|- 4j events at ^/s = 7 TeV and ARjj = 0.7. 

The right hand side of Fig. 1 presents the total visible mass distribution in Z -|- 2j 
production with the same energy and angular separation. It clearly shows that MPI events 
are produced with a smaller center of mass energy than SPI ones. Whether or not these 
different kinematical distribution can be exploited to further increase the MPI fraction 
in the event sample depends on the behaviour of the additional radiation produced in 
association with the hard scattering(s) which is bound to distort both the total visible 
mass and the relative orientation of jet pairs. A dependable estimate of these effects require 
to pass the hard events to a showering Monte Carlo, keeping in mind the normalization 
uncertainties mentioned in Sect. 2. 

The only MPI mechanism contributing at tree level to 7 -|- 3-jets is the production 
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of two jets in one scattering and of a photon and a jet in the other. Therefore, when the 
photon threshold is large, a jet of comparable transverse momentum is also present. This 
feature could reasonably be expected to provide an additional tool to significantly reduce 
the SPI contribution. Unfortunately, as shown in Fig. 2, only a modest reduction can be 
achieved in this way since the pT spectrum of the highest transverse momentum jet is quite 
hard in SPI events. 

4. Conclusions 

In this paper we have estimated the contribution of Multiple Parton Interactions to Z + 
2/3/4-jets and 7 + 3-jets production, comparing the traditional factorized double par- 
ton distribution functions, using both MSTW2008LO and CTEQ6L1 PDF, and the new 
correlated set by Gaunt and Stirling. 

The predictions for the GS09 correlated dPDF differ by at most 15% from those with 
MSTW uncorrelated distribution functions. The uncertainty due to the choice of PDF is 
in the 30 to 90% range. 

It appears quite feasible to measure the MPI contribution to Z-l-2/3/4 jets already 
in the first phase of the LHC with a total luminosity of one inverse femtobarn at 7 TeV. 
If as expected the trigger threshold for single photons is around 80 GeV, the Z + 2-jets 
process may well turn out to be more easily reachable than the 7 + 3-jets channel. It is 
worth recalling that the results presented here are expected to be valid also for + 2/3/4- 
jets with a larger cross section. The possibility of measuring the MPI fraction in several 
channels could allow to extract double parton distribution functions from the data. 
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